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1. INTRODUCTION 

The main purpose of this grant is to investigate the application of 
modern control theoretic ideas to the design of flight control systems 
for the F-8 aircraft. More specifically, the design of an adaptive 
control system based upon the so-called multiple model adaptive control 
(MMAC) method is under consideration. 

In this progress report we shall discuss in an informal \my the 
progress so far. Technical details will not be given since they can be 
found in a sequence of informal interim reports transmitted to NASA 
Langley Research Center as well as an oral presentation to Mr. Elliott 
(the grant supervisor) and his associates at NASA/LRC. In addition a 
partial summary paper is included as Appendix to this report. 

2. SUMMARY OF CONCLUSIONS AND FUTURE PLANS 

A great amount of insight has been gained on the MMAC algorithm 
during this time period from both a theoretical and a pragmatic point 
of view. We summarize these below. 

2.1 Convergence of the MMAC Algorithm 

The transient and steady state behavior of the MMAC algorithm 
depends upon two factors 

(a) the relative sizes of the residuals generated 
by the banks of Kalman filters in real time. 

(b) the prior dynamic and statistical information 
associated with the sensor noise and ^^ind 
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disturbance variances, which is £i\ir:warizod in 
the residual covariance matrix which does not 
depend upon real time measurements. 

The effect of both of the above factors is extremely complex, and 
difficult to understand from an analytical point of view. Nonetheless, 
a rearrangement of the dynamic stochastic equations that govern the 
probability evolutions (see Appendix A) has provided insight into the 
problem. Prior researchers have claimed tliat if the correct model is 
included in the bank of models, then it will be identified with prob- 
ability one. We believe that this statement is only valid under certain 
additional assumptions (e,g, the existence of some persistent excitation) 
and that a rigorous proof in the MMAC context is not available as yet. 

We are in the process of investigating under what precise conditions 
one can guarantee the convergence and stability of the MMAC algcrithm, 
as well as introducing persistent excitation inputs to the aircraft. 

2.2 Sensitivity of the Ml'lAC Algorithm 

The implementation of the MMAC algorithm indicates that its transient 
response can be quite sensitive to both plant and sensor noise. The 
probabilities tend to jump very rapidly, and very little probabilistic 
averaging is evident. This is due to the great ’'amplification” inherent 
in the exponential functions in the equations that dynamically propagate 
the probabilities. This mathematically "optimal" behavior can be 
undesirable for aircraft implementation. For this reason, we are 
currently modifying our MMAC subroutines so that after the probabilities 
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are generated in the "optimal" v/ay, a moving average {whose window*»lChgth 
V7ill have to be determined) of the probabilities will be generated, 
and it is this moving average probability which will v/eigh the controls 
generated by the LQG compensators. V?e conjecture that this simple 
moving average procedure will cause the control system to be less 
sensitive to noise. 

2 . 3 The Value of Lateral Measurements 

In our previous status report we conjectured that the value of 
information provided by the lateral sensors in the overall MMAC prob- 
ability calculation was outweighed by its computational burden. At 
this stage of our research we feel that this conjecture has to be 
reexamined. 

One of the reasons that in our earlier simulations (June 1975) 
v;e felt that the lateral sensors did not provide sufficient information 
for the probability evaluation was due to the fact that the design of 
the lateral Kalman filters included a large amount of "fake" driving 
white noise (to partly compensate for the difference between Gera's 
and Wooley's models). As a consequence/ the lateral Kalman filters 
had very high bandwidths. This caused the selectivity of the lateral 
Krilman filters to decrease when incorporated in the MMAC design. To 
put it another way the lateral MMAC design had difficulty in distinguishing 
models; as a result the probabilities could result in mismatched unstable 
combinations . 

Once this problem was recognized, we redesigned all lateral Kalman 
filters. Our current design (tested at LRC in September 1975) is based 
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upon 


(a) Wooley's data for the linear open-loop 
dynamics 

(b) The only effective source of plant white 
noise is that associated with the wind 
disturbance. 

As a result of this change wo observed from the simulation results 
that the lateral dynamics can provide valuable information with 
respect to model identification, relatively rapidly, and that the 
MMAC design is less sensitive to sensor noise. Effectively, the 
new lateral Kalman filters have a smaller bandwidth, and they are 
more selective. 

In general, when evalup'-'?d in a regulator context the lateral 
MMAC system works very v;ell. 

2.4 Combined Probability Evaluation 

As explained in prior status reports, as well as in the paper 
reproduced in Appendix A, the longitudinal and lateral probabilities 
are generated separately, and if desired they can be easily combined 
in an overall set of probabilities. Once more, based on "theoretical" 
considerations, this combination of probabilities appears to be the 
correct approach. 

There are potential difficulties, however, with this approach 
which require further investigation, due to the fact that the prob- 
ability transitions tend to be very rapid. For the set of flight 








conditions that v/e have available, one can construct two stability 
mismatch tables, one for the longitudinal dynamics and another for 
the lateral dynamics. Effectively, what the miamatch stability 
table indicates is whether or not the coinbination of the j"th LQG 
compensator in feedback about the i-th flight condition is stable. 

The difficulty arises because the stability mismatch tables are 
quite different for the longitudinal and lateral dynamics. In the 
simulations conducted at LRC in September 1975 we observed that quite 
often the MMAC control system "worked" better if controlled from the 
individual probabilities rather than the combined ones. This may be 
due to the rapid changes in the probabilities and the problem may be 
alleviated if the "moving average" scheme discussed in 2,2 is imple- 
mented, We plan to continue to investigate the desirability of using 
the combined probabilities in the overall MMAC design. 

2,5 Pilot Input Design for Longitudinal Dynamics 

During this time period two different designs for handling pilot 
inputs were obtained and tested at LRC. 

The first design translated 

(a) stick position 

(b) stick rate 

into commanded values of the longitudinal state variables using first the 
linearized short period dynamics, and then switching to the long period 
dynamics. This design was tested at LRC in June 1975, It was found 
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unsatisfacbocy because tlie difforonco in steady state values obtained 
from the linearized dynamics are significantly different from those 
obtained from the nonlinear dynamics. 

The second design used exclusively the short period dynamics 
and translated stick position into a commanded pitch rate and normal 
acceleration. These values were then used to generate error signals 
which were fed iivho the existing WIAC regulator control system. This 
design was tested with pilot inputs (provided by C. Wooley) at LRC in 
July 1975. The response was deemed acceptable ^ although somewhat sensi- 
tive (higJi gain) . 

In the short term we plan to still use this later design. However / 
it is not clear how one can evaluate it in terms of the C* criterion. 

For this reason# we are currently carrying out a feasibility study in 
which a C* - like design is obtained, so as to evaluate its impact on 
the overall MMAC concept, 

2,6 Pilot Input Design for Lateral Dynamics 

The approach used for the lateral system was distinctly different 
than that used for the longitudinal system. For the lateral system the 
pilot inputs were introduced to a linear model (provided by LRC) which 
generated the desired responses for the lateral state variables. 

First, we wished to examine if a very simple scheme would work. 
Toward this go?:l we simply used the model generated state variables as 
the desired time varying reference trajectories, and then generated a 
set of error signals which were introduced in the MMAC system. This 
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doBign was toctcd at LRC in Juno. It's responso v/ag judged unsatisfactory 
in view of the largo commanded lateral accelerations. 

At the present tl.ma we are nearing completion of a true model- 
following design which should alleviate the high lateral accelerations 
associated with the prolitni.nary design. 
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SUMMARV OF REAL TIME SIMULATION 


EXPERIMENTS AT LRC 

Trip on 6/4/75 •* G/6/75 

Tasks accomplished Before the Trip ? 

1. Finished tho computations of new regulator and i'ilter gains for 
the lateral dynamics and implemented in the control progrcim. 

(Reason of this change of design has been discussed in our oral 
report on 7/23/75.) 

2. A rudimentary pilot command system for both longitudinal and 
lateral dynamics was designed and implemented in the MMAC control 
system. These designs were carried out under two different 
approaches ; 

For the longitudinal dynamics; A steady state response 
following technique was used. 

For the lateral dynamics; A '’model following" type of 
technique was used. 

3. The procedure of inputlng MMAC system data was changed so that 
the core requirement of implementing the MMAC system was reduced 
significantly. All system data was stored in a tape ( or a fast 
memory device) and only those system matrices used on the real time 
computation for the ^JMAC scheme were called and stored in the pro- 
gram during the reset mode. 
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Simulations Dono; 


1, In order to answer the conjecture posed in the last semi-annual 
report, that the information obtained from lateral dynamics was 
not important, we have done the following experiments for P/C #19 
(1.4 Mach, 40k ft> and F/C #11 (.6 Mach, 20k ft) with different 
combinations of models in the filter banks 

(a) V?ith initial conditions on 2® sideslip: 

1) Both longitudinal and lateral dynamics m-lAC systems 
were operating; 

2) Only lateral dynamics MMAC system was operating; 

3) Probabilities were updated b/ t,;e longitudinal 
information only, and the lateral control was 
computed by the following two different approaches: 

(i) Averaging the "optimum" controls computer 
from the outputs of each filter in the bank. 

t 

(ii) Averaging the controls computed from the 
outputs of a single filter. 

(b) With initial conditions on 6° angle of attack and 2“ sideslip. 
(Repeat the same experiments in (a) . 

(c) With thunderstorm (Repeat (a) and (b) ) . 

2, Similar experiments were carried out for the lateral pilot command 
system. 

3, Single model pilot command system for the longitudinal system was 
tested for F/C #7, F/C #11 a..d P/C #19. 


4 
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Soaio ConcluBlons front thoso Exporlmoiita ; 


1. The new regulator design for the lateral dynamics performed 
well and the xjrior InsbabilitieB on the bank angle response 
never occurred. 

2, There was no significant improvement on the overall system 
performance while the lateral dynamics information was used. 
(It has been found since then that this conclusion was in- 
correct due to excessive fake white noise in the filter 
designs.) 

3, The lateral pilot command system generated a higher lateral 
acceleration than the value generated by the model which was 
due to an ”ad hoc" approach on the model following design. 

4. The longitudinal pilot command system did not v;ork as well as 
we expected. There v:ero some problems on the switching logic 
which controls the pilot coiwand system to follow steady state 
values for the short period response and the long period 
response, rcspoctively. 
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Trip on 7/22/75 - 7/23/75 

Tasks Accompllahed Doforo tho Trip t 

1. To tosb tho tentative conclusions we had reached in the June 
trip, the entire lateral dynamics filters wero rnmovod from the 
control program. Only a single Kalman filter was implo'jiented 
for the lateral dynamics, tlio control was computed from state 
estimates of the filter and tho optimum regulator gains for each 
model in the MMAC bank and then averaged by the probabilities 
computed from longitudinal dynamics. There were 7 models in tho 
filter bank instead of the 4 models used on the June trip. 

2. Tho longitudinal pilot command system was redesigned by eliminating 
the long period response following scheme. 

Simulations Done ; 

1, Complete tost of the now system. 

2. Flight the control system from the cockpit. 

Conclusions ; 

1. Xiongitudinal pilot command system was acceptable, although 
the system was a little bit sensitive to the stick command 
(but it can bo fixed by scaling down the gains on the ADCH 
channel input) . 

2. Lateral pilot command system had instobilities after 20 seconds 
of simulations which happened because of mismatch instability 
and was not observed on the June trip because we did not run 
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the system that long. Vo fix th.'ls problomr wo need a corroctod 
filter design for tlio lateral dynamics and uno both longitudinal 
and lateral dynamics information for the MilAC schema. 

Trip on 9/2/75 - 9/5/75 

Tasks Accomplished Before the Trip t 

1, Chuck Wooley’s linearized model was used for the filter design, 
unnecessary fake jioises for the modeling error v/ere eliminated 
and the wrong turbulence noise statistics was corroctod. 

2. The complr-'t;*? system was implemented. 

Simulations Done and Conclusions ; 

Over three hundred experiments for the MMAC regulator system 
and pilot command systtm for the longitudinal dynamics were done 
during this trip. Several new findings were obtaineJ during this 
trip) we arc still on the stage of analyzing this massive data; 
a detailed report will be written at a later date. 
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Abstract 

Ibo purpose of this paper is to suiranarlno 
results obtained for the adaptive control of the 
F-BC aircraft using the so-called MMAC method. 

Tlie discuBsion includes the selection of the per- 
formance criteria for boOi the lateral and the 
longitudinal dynamics, tho design of the Kalman 
filterB for different flight conditions, tho 
'Mdentifiention'' aspects of tho design using 
hypothesis testing idoos, and tho performance of 
the closed loop odaptive system. 


1. Introduction 

Tho purpose of this paper is to present 
preliminary results on a study which involves 
tho application of advanced adaptive control 
techniques to the design of a stability augmen- 
tation system in both tho longitudinal and 
lateral dynamics of the F-BC aircraft. NASA has 
been using the F-OC aircraft as s test vehicle 
for evaluating different digital-fly-by-wire 
{DFDW) contrc.l tochniques, using the IBM AP-101 
as the airborne computer. VJo remark that tJie 
eventual implementation of the control algorithms 
on tile specific airborne computer hos hod n major 
impact upon tho philosophy adopted for the design 
of tho control system in view of the obvious sto- 
rogc and real-time computational constraints. In 
addition, the design was crucially dependent upon 
the sensors that could be utilized in the sense 
that sensors that utilized external aerodynamic 
moesuremonta , e.g., airspeed, altitude, angle of 
attack and sideslip vanes should not be employed 
in the cond.date design. Thus, the design guide- 
lines required that the sensors associated with 
the adaptive control system should be limited to 
accclcromnkers, rate gyros, and perhaps attitude 
sensors (.ilthough the latter were deemed undesir- 
able in view of their errors when the oircraft 
underwent aeverc pilot induced maneuvers) . 


•The theory ond initial olgorithm deveXopment 
associated with this study were developed with 
support from NASA/Ames Research Center under 
grant NGL- 22-009-1 24 and from AFOSR under grant 
72-2273. The specific application to the F-OC 
was supported by NASA Langley Research Center 
under grant NSC-IOIQ. 


From tho viewpoint of modelling, it is ob- 
vious that tho dynamic state equations of an air- 
craft involve nonlinear differential equations 
(see Etkln (11). However, the information given 
by NASA Langley Research Canter (LRC) to tho MIT/ 
ESL team consisted in tlio specification of the 
uncoupled , linear time-invariant open-loop 
longitudinal and lateral dynamics of the F-8C 
aircraft associated with equilibrium flight . 

Table I gives a list of the flight conditions 
that were available for the design. Titus, the 
general structure of the equations wore of the 
form x(t)»>^(t)+Du_(t) . Tho numerical values 
of tho elements of tho A ond matrices can be 
found in a report by Gera (2), based upon wind 
tunnel tests, and o report by Wooley and Evans 
[31 , baaed upon linearization of the nonlinear 
dynamics employed by NASA/LRC for tlteir nonlinear 
simulation of tho F-BC aircraft. We remark at 
this point that tho numcricoi values for the j\ 
and B matrices given in (2l and (31 are not iden- 
tical reflecting the fact that different sources 
wore used to obtain them. The design reported in 
'this paper is based upon Gera's report (2]. 

The fact that the 16 flight conditions span 
an extremely wide envelope for operating tho 
aircraft, with drastic changes in the open-loop 
dynamics, makes tho fixed-gain design of tho con- 
trol system unrealistic. Furthermore, handling 
qualities requirements, ouch as the C* criterion, 
Indicate that pilots desire different closed-loop 
dynamics at different flight conditions. Thus, 
some sort of "adoptive” gain-scheduling control 
system was required. However, straight-forward 
gain scheduling based upon quantities such as 
velocity, altitude, and dynamic pressure was not 
permitted in View of the sensor restrictions men- 
tioned above. Hence, the adaptive control system 
had to be designed in a novel way. 

An additional restriction on the design was 
that the sensor noise and wind disturbances had 
to be incorporated. This led to the need for 
employing Kalman filters, with constant coeffi- 
cients because of tho computer memory limitations . 

The above problem overview sets the ground 
for the specific adaptive control technique which 
we selected to investigate in great detail, vre 
call the adaptive control technique the Multiolc- 
Hodcl-Adaptivo-Control (MMAC) method, and wo shall 
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(Uscups it in ntoto dotail in Soction 5 of thlt) 
pnpor. It 1b only ono of nayoral tothnlquoo 
bnood upon dovoloprnunto in tnodorn control theory 
(doc the survey article by Athana and Varaiya 
(4]) and it hao ita origins in combining 
hypothoBlD-tooting and otochaatle control ideas 
(boo roforonccs 141 to (Q])< It was eolectod for 
thlu ctudy because of ita potential promiao in 
academic cxaitiplca tSI-18], and because its memory 
and real-time conipubatio.'al requirements could bo 
readily assessed in View of its non-ltorativo 
nature • 

An explained in more detail in Soction 3 , 
the UMAC method requires that a full blown n toady 
ntato Unoar-CUadrotic-Gausaian (tOG) controller 
bo implemented for each flight condition. This 
necoBDltatcd the development of suitable quad- 
ratic performance criteria for both the longitud- 
inal and the lateral dynamics) those are 
doocrlbod for the continuouo time case (91 in 
Sections 2 and 3, respectively , For implementa- 
tion, one needs a discrete time tQC controller 
(lOl, This is described in Section 4, together 
with the discussion of sensor errora, The HWAC 
algorithm is described in Soction 5. The simula- 
tion results using the nonlinear F-OC dynamics 
are described in Soction 6, Soction 7 presents 
the major conclusion of our studios so far. 

Wo remark that in this paper wo shall only 
focus our attention to the regulotion aspects of 
the problem, i.e,, return to equilibrium flight 
from some Initial conditions and in the pre- 
sence of stochastic wind disturbances, In our 
study we arc considering the proper way of incor- 
porating human pilot Inputs for both the longi- 
tudinal and lateral case. However, we shall not 
present in this paper any of the approaches and 
preliminary results for the pilot input case. 


2. Longitudinal Dynamics 

2.1 Introduction 

In this section we present an overview of 
the LQG philosophy adopted for designing the reg- 
ulator for the longitudinal dynamics. Attention 
is given in the development of the quadratic per- 
formance index and the subsequent model simplifi- 
cation using a short period approximation. The 
main concept that wc wish to stress is that the 
quadratic porformancu criteria employed changed 
in a natural way with each flight condition. The 
surprising result was that the short period poles 
of the resultant longitudinol closod-loop system 
v/cre characterized for all flight conditions by 
two constant damping ratios, one associated with 
all subsonic flight conditions and one associated 
with all supersonic flight conditions. 

2.2 The Longitudinal State Description 

Because of a rate constraint saturation on 
the elevator rate, the control variable selected 
was the time rate of change of the commanded 
elevator rate (fQ,.(t)). This was integrated to 
generate the actual commanded elevator position 


which was introduced to a first order 
servo with a time constant of 1/12 seconds to 
generate the actual deviation of the elevator 
6g(t) from its trlmned value. The elevator was 
then related to the four "netural’* longitudinal 
state variables namely pitch rate, q(t) (rad/soc) , 
velocity error v(tMft/Bec>, perturbed angle of 
attack from its trimmed value, a(t)(rad), and 
pitch attitude deviation from its trimmed value 
6(t)(rad), (2), In addition, a wind disturbance 
state w(t) WAS included (see Appendix A), Thus 
the state vector x(t) for the longitudinal dynam- 
ics was characterrzod by seven components 

(2.1) x'Ctl-tqCt) , v(t), Ci(t), 0(t), i^(t), 
«„„(t), w(t)J 

OC 

and the control variable u(t) was the corttnanded 
elevator rate 

(2.2) u{t)-4^^(t) 

This led to a linear-time invariant characteriza- 
tion for each flight condition of tho form 


(2.3) x{t)>»A^x(t)+nu(t)+I^C(t) 

where ^(t) was zero mean whlto noise, generating 
the wind disturbance and accounting for random 
actuator errors. The elements of ^ and 
chonged with each flight condition while 

(2.4) B-(0 0 0 0 0 I 0)' 

2.3 Tho Longitudinal cost Functional 


In order to apply the standard steady state 
LQG procedure (9) a quadratic performance index 
has to be selected. The general structure of the 
index was 


(2.S) a-/^'(t)2iX(t>+u'(t)R^U(t)dt 

Note that the weighting matrices g, , had to be 
different from flight condition to'flight condi- 
tion reflecting in a natural way that the pilot 
wants different handling qualities as the speed 
(and dynamic pressure) changes. 


In the initial design it was decided that 
one should rtO.atc the maximum deviations of 
o pitch attitude, Omax 
o pitch rate, qmax 
o normal acceleration, * 

o maximum connanded elevator rate, £g^max 
resulting in tho following structure of 
the performance criterion . 

„ T*"®**"’ 


nzmax 


dt 


ecmax 


The normal acceleration g's, was not 

used as a state variable. However, it is linearly 
related to some of the longitudinal state vari- 
ables according to the formula 

(2.7) a^^(t)-~jk, v(t)+Kja(t)+K36^(t)j in g's 
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Vt boin 9 tha oqull,ibrium spaod. Tho con&tanto 
ki( k 2 > k) can bo calculated Cirom the open loop 
*■ matricoB, and honco change with flight condl- 
rloni Effectively the etructurc of the criterion 
(2.G) implies that if at t«0 the inaximiin values 
of accoloratlon> pitch cate< or pitch attitude 
occurred^ then one would be willing to saturate 
the elevator rate to remove them, For the pre- 
liminary design the following numerical values 
wore selected (witli the help of T, Elliott and 
J. Cera) 

( 2 . 0 ) 

,6 "0.435 rad/soc 

where aj> is the (3.3) element of the open loop 
longitudinal ^ matri'.<. 

Roughly speahing. this criterion neons that 
one is willing to soturoto tho elevator rote 
(0.435 rod/sec for tho F-GC) if a normal accel- 
eration of Gg was felt, or a pitch rate equiva- 
lent to log's, or a pitch error which if 
translated to angle of attock would also generate 
a Gg normal acceleration. 

Tl)c abovQ numerical values were translated 
into tho appropriate matrix (non diagonal pos- 
itive Gomldcfinitc) which changed from flight 
condition to flight condition, while R,“R;*1/ 
(0.435)^ Cor all flight conditions. Hence, the 
resulting *J?G problem could be solved using 
available computer subroutines [11]. 

2.4 Reduced Longitudinal Design 

Tho design was modified for two reasons. 

First the gain from the velocity state varioble 
v(t) was extremely small. Second, it was desir- 
able to avoid using the pitch sensor. The pitch 
0(t) is weakly observable from the system dynam- 
ics BO that even if a Kalman filter was used in 
the absence of pitch mcasui'cments, large estima- 
tion errors would be obtained which would adversly 
affect tho performance of the control system since 
there is significant feedback from the estimated 
pitch attitude. At any rate, since a pilot would 
fly the aircraft he would be able to control 
pitch himself. 


tive tradeoff between tho maximum normal accelera- 
tion, and maximum pitch rate, This 

is consistent with tho C* criterion (12). 

When tho eliort-period closed-loop poles wore 
evaluated for both designs using the numerical 
values given by (2.B^, wo found tho unoxneeted 
result that the damping ratio was constant (0.4B0) 
for all 11 subsonic flight conditions, and also 
constant (0.3G1) for all the supersonic fllglit 
conditions , Tho closed-loop nafurolly frequency 
Incroasod with dynamic pressure « 

Since no poio-placcmont techniques wore em- 
ployed (l.e., the mathematics were not told to 
place tho closed-loop polos on a constant damp- 
ing ratio lino) , we constructed a tradeoff by 
changing (decreasing) tho maximum pitch rate 
“Imax ' would increase the pitch rate penalty 

in tho cost functional, and one would expect a 
higher damping ratio. The following values of 
«lmax “ote employed 

(2.10) , DgAo, 6g/V(, 4g/Vo 

Once more tho constant damping ratio phenomenon 
was observed, i.o, , for each value of q^^^ tlic 
short period closed loop poles for all subsonic 
flight conditions foil on a constant damping ratio 
lino, and similarly for all supersonic flight 
conditions. This was further verified by consid- 
ering an additional 13 different flight condi- 
tions. 

The numerical results are presented in Table 
II, The reason for this regularity of the solu- 
tion of the LQ problem is under Investigation. 


3, Lateral Dynamics 
3,1 Introduction 

In this section we present tho parallel 
pSiilosophy for the development of tho control 
system for the lateral dynamics. In this case 
the development of o p<»rfomanco criterion was 
not as straight-forward as in the case of the 
longitudinal dynamics. For an extensive discus- 
sion sec the s.M. thesis by Greene [13], 


This led us to eliminating the velocity error 
v(t) and pitch 0(t) from the state equations and 
obtaining tho "short period" approximation (5 
state variables). Since pitch did not appear the 
criteria (2.6) was modified to 
*(t) 


(2.9) d! 


,L0N 




hzmax 


(t) 


q^max J 


GO 

+ ^ 


dt 


ecmax 


and the resultant LBG problem was resolved. 


2.4 Summary of Results 

From the viewpoint of transient responses to 
the variables of interest (normal acceleration, 
pitch r.ato, angle of attack) tho transient re- 
sponses to initial conditions were almost identi- 
cal for both designs. Thus, the short period 
motion of tho aircraft was dominated by the rola- 


3,2 The Lateral Dynamics state Model 

The control variables selected for lateral 
control were 

(3.1) u I (t)"^_(t) "commanded aileron rate 

* (rad/scc) 

(3.2) ui(t)«i (t)'conmandod rudder rate 

(rad/sec) 

so that the control vector is defined to be 

(3.3) u’{t)-[ui(t' ui(t)l 

Tho servomechanici wore taken into account. Tho 
conmanded ailero'i and rudder rates were integrat- 
ed to generate the commanded aileron (^,^c(^) ) and 
rudder (6j.g(t)’ positions, respectively. For the 
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r-OC aircroCt the commanded alloron rate 
drtvoa a first order la^ oorvo, with a time con- 
struct of 1/30 aecondOf to generate the nctuoX 
aileron iwoltlon (rndo) . 

Tlio commanded rudder rate (rads) drives 

a first order lag sorvo< with a time conntant of 
1/25 seconds* to generate the actual rudder posi- 
tion drCcHrads). The actual aileron and rudder 
position* dci(t) and dj-{t) * then excite the four 
"natural" lateral dynamicH state variables* name- 
ly roll-rate p(tMred/sec} * yaw-rate r(t)(rad/ 
sec)* sidoollp angle {)(t)(rad), ond bank angle 
(^{t}(rad). In addition* o wind disturbance state 
variable w(t)* soe Appendix A* drivoe <tie equa- 
tions in the same way as the sidosllp variable. 

lliua* the state equations for the lateral 
dynamics are charncterired by a 0-dimcnsional 
state vector x(t) with components 

<3. 4) x'(t)-[p{t) r(t) P{t) d(t) djit) d^(t> 



acmax rcsuuc 


The following maximum values were used 

Maximum lateral acceleration, aynyix^O'^Sg's 

4Vd 

H*Klmum roll rAte^ p (a*..kio) 

"“’‘v/Iog 

Maximum sideslip angle* 0 a-_ 

'"“v/Io'g 

Maximum bank angle, (“45*1 

Maximum coitinandod aileron rate»l.C3 rad/sec 
Maximum coneianded rudder ratc<*l«22 rad/sec 


d„^(t) d^^tt) w(tn 

ac rc 

and the overall lateral dynamics toke the form 


(3.5) K<t)"Aj_x(t>+B u(t)+b^^(t) 


whore the rcro mean whitu noise vector ^(t) gen- 
erates the wind disturbance and compensates Cor 
modelling errors. Once more the matrices ^* 
change with flight conditions (2) , (13) while 


(3.6) 


> fo 00000100 
-"[o 00000010 


3.3 The Lateral Cost Functional 


ibc lateral performance index used (after 
Boverai iterations) weighted the following 
vurioblcs I 


Gee (13] for an extensive discussion of how this 
performance criterion was dorivod; and B33 
are obtained from the open loop ^ matrices 

Tliore is no natural way of arriving at a 
simplified model for the lateral dynamics* as was 
the case with the longitudinal dynamics. Hence 
the bank angle cannot bo eliminated. Although a 
bank angle sensor was deemed undesirable , the 
weak observability of the bank angle caused largo 
state estimation errors* using Kalman filters* in 
the bank angle and the sideslip angle if a bank 
angle sensor was 1. t included. For these reasons* 
it was decided to employ a bank angle sensor and 
to penalize bank angle deviation* because bank 
angles larger than 20* can introduce significant 
nonlinearities through trigonometric functions 

(ll. 


o lateral acceleration* a (t) 
o roll rate* p(t) ^ 

o sideslip angle* 0(t) 
o bank angle* ^(fc) 

VG 


(in g's) 

(in rad/scc) 
(in rad) 

(in rad) 


o commanded aileron rate, ^gg(t) 
o commanded rudder rate, ^j;c(t) 


Tlic lateral acceleration* a (t) , is not a state 
variable. However, for small perturbations from 
equilibrium flight, it can be expressed as a 
linear combination of the lateral state variables 
and the trim angle of attack, Qo* by the following 
relation 

(3.7) ay(t)=~ ((kj-ao)p(t) + (k2+l)r(t)+ks0(t) 

+ki,i}.(t)4k}6j.(t))-4.(t) 

where the constants ki*...* ks can be found from 
the lateral open loop ^ matrix, and change with 
the flight condition. 


The following structure of the quadratic 
performance criterion was established, 


Once more* the JJQ can be solved. Notice that 
the use of the pcrfciDnanco criterion (3.8) results 
in a state weighting matrix g, (non-diagonal) 
which changes with flight condition. 

3.4 Summary of Hesults 

The above performance criterion gave reason- 
able responses for a variety of initial conditions. 
Its main characteristic is to reduce any lateral 
accoloratlons (by forcing the aircraft to go In 
coordinated turns) and to null out bank angle 
errors in a slower tronner. 

Once more we observed a constant damping 
ratio ( .515 ) for all supersonic conditions and 

a relatively constant damping ratio { .025 ) for 

all subsonic flight conditions. Ho additional 
tradeoff studies were conducted by changing the 
weights In the cost functional. 
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4i scntsocs, Kalinnn rllteza and 
Dloct^otc M3G Conponontora 

4il Introduction 

Itto digital implomentation of tho control 
system requires the dlacrato-titno solution of tho 
liQG problem (10). As wo shall see in tho next 
ooctionf tho MHAC approach requires tho construc- 
tion of a bank of uyj controllers, each of which 
contains a diccroto Kalman filter (whoso resid- 
uals aro used in probability calculations and 
whose state oBtimatcB aro used to gonorato tho 
adaptive control signals), lienee, in this section 
wo present on overview of tho issues involved in 
the design of the LQG controllers based upon tho 
noisy sensor moosuroments . 

4.2 The Sampling Interval 

A sampling rate of 8 roeasuroments/socond 
was established, such a slow sampling rate was 
selected so os to bo able to carry out in real 
time tho multitude of real time operations re- 
quired by tho Ml’lAC method. 

4.3 Sensors and Poise Characteristics 

As explained in tho introduction, tho guide- 
lines for design excluded tho use of air data 
sensors. Thus, measurements of altitude, speed, 
angle of attack, avid sideslip angle *fere not 
available. After some preliminary investigations 
It was decided that sensors that depend on trim 
variables (elevator angle and pitch attitude) 
should not be used so as to avoid cBtimatlr.g trim 
parameters. Table III lists the sensors and their 
accuracy characteristics that were used in this 
study. We stress that the sonsoro measure the 
true variables every 1/B seconds in the presence 
of discrotc-sero mean white noise with the stan- 
dard deviations given in Table III. 

Finally, we remark tliat in this study wo 
assumed th-*’ all sensors wore located at the C.G. 
of the aircraft. 

4.4 Tho Design of Kalman Filters 

For each flight condition tho steady-state 
discrete-time Kalman filter, with constant gains 
was calculated, for both t)fc longitudinal 
and lateral dynamic models. The level of the 
plant white noise associated with the wind dis- 
turbance generation was selected so that wo 
assumed that tho aircraft was flying in cumulus 
clouds. (Sec Appendix A.) 

The decision to use steady state constant 
gain Kalman filters was made so as to minimize 
the computer memory requirements. 

Finally, we remark that in view of tho slow 
sampling rate, the continuous time filtering' 
problem was carefully translated into tho equiva- 
lent discrete problem ( 13) to (15). 

Tho constant covariance matrices of the 
Kalman filter residuals, denoted by §lLAT 


for the longitudinal models end lateral models 
were computed for each flight condition denoted 
by 1. Ae we shall see these are Important in tho 
generation of the MHAC variables. 

4.5 The Design of tho Discrete 14^3 Compensators 

Through the use of the separation theorem one 
con design the discrete DQG compensators. This 
implied that tho 14 ? problem defined in continuous 
time in Sections 2 and 3 had to bo correctly 
transformed into the equivalent discroto-tlme 
problem in view of the slow measurement rate, 
Effectively, wo have used the transformations 
given in roforonceo (11), (13) to (IS). 

4.6 Recapitulation 

For each flight condition, indexed by i, a 
complete discrete-time, steady state, 14 ?C compen- 
sator was designed for both the longitudinal and 
lateral dynamics. Each componoator generated 
every 1/B second tho optimal control, namely the 
optimal eomandod elevator rate 6^. (t) for the 
longitudinal dynamics, and the optimal oammanded 
aileron rate o,|e(t) and rudder rate ^j.(.{t) , based 
upon the noisy measurements of tho appri;,priatc 
sensors (See Table lU) every 1/S second. 

Docause of tho appropriate transfc'rmations 
of the continuous time IQG problem to the discrete 
one, we noted no significant degradation in per- 
formance at this low campling rate. 

The need for adaptive control is obvious be- 
cause if wo assume that tho aircraft is in flight 
condition i, but we use tho compensator ob- 
tained for flight condition j for feedback con- 
trol, this mismatching may generate either an 
unstable system or, often, a system with degraded 
performance. 


5. The HMAC Method 

5.1 Introduction 

In this section wo present tho basic Idea 
behind the MHAC method, and discuss how it was 
used in tho F-6C context. In particular, we 
demonstrate how tho Information generated by the 
lateral and longitudinal sensors is blended to- 
gether. Finally wo make some remarks associated 
with tho MHAC method and its general applicability 
to the design of adaptive control systems. 

5.2 The Sasic Idea 

Suppose one has H linear, dlscretc-tlmc 
stochastic time-invariant dynamic systems , ^ ndexed 
by i»l, 2, N, generating discrete-time 
measurements corrupted by white noise 
Suppose that at t"0 **naturo" selects one of these 
systems end places it inside a "black box." The 
true system generates a discrete set of measure- 
ments £{t) , The objective is to apply a control 
signal u (t) to the true model. 
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Tlio vorelon of U)o HiVVC motJiod oniployod 1 b 
• fl follows ( one conotructs a dloceoto>'timG stoady 
state 1/QG controller for onch iradolt thUBi one 
hae a bank of N MJG coinponaatorB. /\n shown In 
riguro 1, each LQC componaator In drlvon by tho 
actual control applied to the ayatotitf ^(t ) , and 
driven by the actual noisy mdosuramcnt vuctor^ 
z(t) . Hiero arc two signals of Intorost that 
cTach ><QG componoator generates at time t 

(1) tho control vector ^ (t ) , which would bo 
tho optimal control If indeed the system 
in tho black box (viz. aircraft) waa 
identical to tho i-th model 

(2) the residual or innovations vector 
r£(t) gonoratod by each Kalman filter 
(which is inside tho i-th compensator) 

It turns out that (see roforoncos t'l] . (5)r 
|V)i (d) examplo) that from tho residuals 
of tho Kalman filters one can recursively generate 
N discrete time soquoncos denoted hy P, (t) , 1*1, 

2 N, t»0, 1, 2, ,..i which under suitable 
assumptions arc the conditional probobilitlos at 
time t, given the past measuromonts ^(T) . x<t 
and controls u_(o) . o^t-l, that the i-th model is 
the true one. 

Assuming then that these probabilities arc 
gonoratod on-line (the formula will bo given 
later) and given that each tQG compensator gen- 
oraton tho control vector uj^(t), then os shown in 
Figure 1> the F!MAC method computes tho adaptive 
control vector u(t), which drives the true system 
(viz. airorn£t)“and each of tho Kalman filters 
inside the LQG compensators, by probabilistically 
weighting the controls ^ (t) by tho associated 
probabilitloB, i.e., 

N 

(5.1) u(t)-j: p, (t)u,(t) 

i-1 ^ 

5.3 Calculation of tho Probabilities (t) 

We assume that at t»0, i.e., before any 
measurements arc obtained, one has a sot of prior 
probabilities 

N 

(5.2) Pi (0), ..., P„(0), P, (0)>0, J:p.(0)-1 

" ^ 1-1 ^ 

that represent our "best guess" of which model is 
indeed tho true one. 

in our version of the MMAC method we have 
available the etoady-Btatc (constant) covariance 
matrix ^ of the residuals associated with the 
1-th Kalman filter. These N residual covariance 
matrices are procomputablo. I«'c r denote the 
number of sensors; then we can precompute the 
N scalars 

(5.3) 0^»«[(2Tt)'^dot 

From the residual vector r^(t) generated by each 
Kalman filter wo generate on-line tho N scalars 

(5.4) m^(t)-r^'(t)s^-*r^(t} 


Then the probabilities at tine t, P^(t), i-1, 2, 

N are comiiutod rocuroivoly from the probabil- 
ities at time t-1, Pj(t-l), by the formula 

p, (t-l)l •oxpC-ffl, (t)/ 2 ) 

(5.5) P^{t)“i;i 

E P. (t-l)fl/0Kp{-ii. (t)/2) 
j-l J 3 3 

with the initial probabilities, P.(0) given. It 
has been claimed that (5), (G) , (0), under suit- 
able assumptions that asymptotically the true 
model is Identified with probability 1. 

5.4 Important Pemarks 

1) It has been shown by Willnor (G) , that 
tho KMAC method, i.e., generating the control via 
(5.1) is no;^ optimal \lt is optimal under suitable 
assumptions for the last stage of the dynamic 
progronmiing algorithm) . 

2) Tho KMAC algorithm is appealing in on ad- 
hoc way because of its fixed structure and because 
its real-time end memory requirements are readily 
computable. 

3) In the version used in thin study, because 
wo use steady-state Kalman filters, rather than 
tine-varying Kalman filters, tho 

P,(t) are not exactly tho conditional 
probabillt'as. 

4) We have been unable to find In tho cited 
litorat!\te a rigorous proof of convergence of the 
claim that indeed tho probability associated with 
tho true model will asymptotically converge to 
unity. 

5) From a heuristic point of view, the re- 
cursive probability formula (5.5) mokes sense with 
respect to identification. If the system is sub- 
ject to some sort of persistent oxitation, then 
one would expect that the residuals of the Kalman 
filter associated with tho correct modal, say the 
i-th one will be "small," while tho residuals of 
the mismatched Kalman filters (j^i, j-l, 2, .... 

N) will be "largo." Thus, if t Indexes tho 
correct model we would expect 

(5.6) n^(t)«»j(t) all j/i 

If such B condition preslsts over several measure- 
ments, the analysis of (5.5) shews that the 
"correct" probability P£(t) will increase while 
the "mismatchijd model" probabilities will decrease. 
To see this one can rewrite the formula (5,5) as 
follows 

FN 

(5.7) Pj^(t)-P^Ct-l)-| E Pj (t- 
Pjl^(t-l) |(1-P^ (t-1) ) &»oxp{-m^ (t)/2)j 
- E P it-l)P;cxp{-m. (t)/2) 

if I ^ ^ 

Under our assiaq>tions 


P»exp(-mj (t)/2) 


il-' 
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(5.0) oxp(-m^(t)/2)ssi. 

(5.9) OKp(-iHj(t)/2)s50 

llonco tho correct probability will ^row accoccitng 
to 

P. (t-1) (l-l*. (t-llp.* 

(5.10) P^(t)-P^(t-l) a. - ^ >0 

t P . ( t-1 ) B. •0Xp{-Bl . ( t) /2 } 
j«l 3 3 3 

which domotiBtratos that as Pj,(t)'^l( the rato ot 
growth olowe down. 

On tho othor hand, for tho Ineorroct modola, 
Indexed by i/L, tho sane assumptlona yield 

-pjt-l)p, (t-l)O.* 

(5.11) P.(t)-P.(t-l)S -3 i i <0 

3 3 

Z P^(t-l)0»oxp(-iaj^(t)/2} 


so that tho probobilitloa docroaoo. 

The same conclusions hold if uc rewrite 
(5.7) in tho form 

(5.12) P^(t)~P^(t-l)-|^ T. Pj(t-l)B*oxp{-mj(t)/2>j''l 

[p. (t-1) Z P. (t-l)/3,*cxp{-m. (t)/2) 

Li 3 1 

-0^*oxpj-mj(t)/2{j 


The above discueeion points out that this "identi- 
fication" scheme Is crucially dependent upon tho 
regularity of the residual boliavlor between tho 
"matched" and "mismatched" Kalman filters. 

(6) The "identification" scheme, in terms of 
the dynamic evolution of tho residuals will not 
work very well if for whatever reason (including 
errors in the selection of the noise statistics) 
tho residuals of the Kalman filters do not have tho 
above regularity assuinptlons. Tobe specific, sup- 
pose that for a prolonged sequence of measurements 
the Kalman filter residuals turn out to be such 
that 

(5.13) mi (t)‘>m 2 (t)«. . .*>m„(t) 

N 


Then 

(5.14) exp(- roj^(t)/2)“ “ for all i 

Under these conditions and (5.12), wc can seo 


that 

1 

(5,15) pj, (t>-P^{t-l)« 


(t-1) Z P.(t-l) (0 "-B *)a 
i/i ^ ^ ^ 

N 

Z PMt-l)0 *0 

j-1 ^ i 


P.(t-l) Z (0.»-0/)PJt-l) 

J J 

N 

Z P.(t-1)0.* 

M J ^ 

suppose that it turns out that one of the 0,*'s, 
and to bo spocifle 0j^*, la donlnant, l.e,, 

(5.16) B,j*>0i* all IrOc 

In this case, tho RiiS of eq. (5.15) will be nega- 
tivo for all iA, whidi means that all the P. (t) 
will decrease while tho probability P)^ (associated 
with tho dominant 0)^*) will increase. This be- 
havior is very Important, cxpeciallv when one lies 
to tho mathcnaticB , and it has not been discussed" 
previously in tho literature to the best of our 
knowledge, 

5,5 Application to tho F-ec 

Tho KKac method can be used in a straight 
forward manner using either the longitudinal or 
lateral dynamics of the F-OC aircraft since we 
have designed both longitudinal and lateral bCG 
compensators for tho available flight conditions, 
as wo remarked in Section 4. 


On the othor hand, wo obtain independent 
information from the longitudinal and lateral 
systems for the same flight condition (i.e. , model) 
indexed by i, Itcnce, it should bo possible to 
blend this combined information into a set of sin- 
• gle probabilities. 

Under the assumption that tho langitudinal 
and lateral dynamics are decoupled K-P Dunn de- 
rived the following relation. 


Let S£ and denote the residual 
covariance matrices of the Kalman filters, for 
the i-th flight condition, associated with the 
longitudinal and lateral dynamics respectively. 
Define 


ei * L 0 N - 
< = ‘“5 ei*LAT- 


^LAT„ 
(2») det£^ 


- 1/2 

- 1/2 


where and are the number of longitudinal 
and lateral senses. Let £l 

denote the Kalaan filter residMl vectors at time 
t, for flight condition i, associated with the 
longitudinal and lateral dynamics respectively. 
Define 


(5.19) 

(5.20) 


“l LOh“£i iio45I*LOM^ LON^*^^ 
*’iiat“^ LAT**^^^ LAT^ IAT* 


Then the overall probability that the aircraft is 
in flight condition i at time t, is generated by 
the recursive formula 


i 



<5.21) 


contributo aono undaratandln 9 upon the khaC nethol 
M ■ doaign .concept. 


The 0* 
to the 

(5.22) 


o xp{ ( t ) /2 ) oxp (-m^ ( t) /2 1|^ 

oxp{-i#jj^jj (t) /2)cxp{-nj jj^^(t) /2)J 

dominance effect dlacuDBod above now rofrra 
relative magnitude of 

»i*-ei*b0/i*WT 


Obvioualy the method ahuuld he expected to worK 
well w)ion both longitudinal and lateral Kalman 
filters arc correctly designed so that the rceid- 
uale of the "matched" Kalman filters ore amallec 
than those of the "mismatched" ones. 


C. Simulation Acsulta 
C.l Introduction 

h variety of almulationo have been done csing 
both a linear model and nonlinear model of the 
F-BC aircraft. Those simulations results arc ty- 
pical. They are solocted such that they con dem- 
onstrate 

1) the speed of identification of the 
MMAC algorithffit 

2) the overall perfomanco of the MKAC 
syotomt and 

3) the 3* dominant behavior discussed in 
Section 5, 


5.6 Discussion 

It should bo immediately obvious, that if the 
MMAC method is applied for the control of the F-OC 
aircraft (oi any other physical system for that 
matter) , one violates a multitude of theoretical 
assumptions. The effect of those upon the perfoim- 
anco of the overall system is difficult to estab- 
lish on an analytical basis, because the MMAC 
system, in spite of its simple structure, ropro- 
Bonts an extremely nonllnoor system. Hence* one 
has to rely on extensive simulation results in 
order to be able to ma ){0 a judgment of the per- 
formances of the overall algorithm. 

Since the alrcra.'t never coincides with the 
mathematical models (‘.ecall the discussion on the 
differences in the data given in rcforancos (2) 
and [31, the P«{t) are not truly posterior proba- 
bilities. Pathor they should be interpreted ns 
time sequences that have a reasonable physical 
interpretation. Honco, in our opinion, the eval- 
uation of the MMAC method solely by the detailed 
dynamic evolution of the Pi(t) is wrong. Rather 
it should bo judged by the overall performance of 
the control system. In the case of the regulator, 
this is easy since one can always compare the 
response of the MMAC system with thst which was 
designed explicitly for that flight condition and 
compere the results, 

Wo remark that such a comparison is much more 
complex vihen one attacks the case of pilot inputs 
which result i,'\ several commanded maneuvers. 

These aspects ate still under investigation. 

There are several unresolved problems as vot 
which pertain to the total number of models to be 
used at each instant of time, how these models are 
to bo selected, how they should be scheduled in 
the absonco of any air data, and how one can arrive 
at a final design that moots the speed-memory lim- 
if-“tions of the IBM AP-101 computer which is used 
in the IIASA F-Bc DFBW program. 

Wo hope that some of the simulation results 
and discussion presented in the sequel can 


Some remarks about the MMAC method are given in 
the conclusions. 

C.Z The Simulation Results 

The simulations wore conducted at a high al- 
titude (40,000 ft), supersonic (Mach 1.4) flight 
condition (F/C *19 in Table I) , Mo plant noise 
was introduced. All models available In the MMAC 
controller wore given equal a priori probabilities 
being the true model. 

Experiment *1 > 

This is a sot of linear simulations with two 
degree sideslip angle ( a B'gust) at time t«0. 

No sensor noise was introduced end the Kalman 
filters wore set at the correct initial conditions. 

Figure 2 shows the probability changes while 
the set of models available in the MMAC controller 
wore F/C 8, 14, IB, 19 and 20. Note that the 
true flight condition was included in the con- 
troller. The correct model is initially chosen 
with high probability within a very short period 
of time (less than 1 sec.) and than switches to 
another model slowly after a few seconds. Lateral 
acceleration is removed within about one second, 
while roll rate and sideslip angle arc reduced to 
aero almost as fast. With no noise perturbing the 
system,. ^ s states of the system have settled to 
near zero afLi:;^ shout five seconds. Thus the 
residuals in all the mismatch stable filters 
approach zero. In this case the 3* dominant be- 
havior discussed in Section 5 occurred. Figure 3 
shows the probability changes when the true model 
(F/C *19) was not included (which was substituted 
by F/C *17) . He observed the some 3 * dominant 
behavior after about five seconds. The most impor- 
tant point to note is that responses of the HMAC 
system arc almost identical. Figure 4 shows the 
responses of lateral acceleration with and without 
F/C #19 in the controller, respectively. 

Similar results wore obtained with other ini- 
tial conditions. However, the speed with which 
the 3 * starts to dominate varies greatly. For 
example, with a roll rate initial condition, it 



haiH’onod much soonor, ilowavori thota la vary 
llttlo dugradation in tiio ovarall ayatgin perfoc** 
Mnco. 

Eaporlwont >2 > 

Thiu la a sot of nonlinaoc almulatlonB with 
an initial olx dcgroo anglo of attack (an 0-guat) . 
Sensor nolao was introduced and the Kalman Clltora 
wore cot at toro initial condltlons< 

Figure 5 shows the change of probobllltlos 
when tliQ sot of models avalloblo in tho KHAC con” 
trollor woco F/C 14, 17, 19, and 20 and when tho 
true flight condition (F/C 19) was included. Tho 
probnbilltioB aro more active than tJioso wo havo 
seen in I'xporimont ttl. It ic bolievod that thn 
fast variation of those probahllltios la duo to 
a combination of tho transiont rosponso of tho 
system and tho noise soguoncoa on tho sonsora. 
Howevor, tho true fllgtit condition io identifiod 
in about 1 aocond. Tho onglo of attack returns 
to its trimmed value within about 2 sccondo, while 
pitch rata and normal accolorotlon ere roducod to 
zero alSiOst as fast. In this case the 0* dominant 
behavior only ocouro Cor a very short period of 
time. Decauao of tho sonsor noise, it is not 
certain if tho drifting of probabllltios aro 
mainly duo to the dominant 8* . Figure 6 shows 
tho probability chsna’on when F/C H9 (true) was 
subBtitutod by F/C Hlli in tho KMAC controller. 
Again, tho rosponscs at tho t!MAC system are 
alirost identical. Figure 7 shows tho responsoo 
of tho angle of attack with and without F/C tll9 
in tho controller, rospoctivoly . 

7. Conclusions 

Based upon many simulations using both a 
linear modol and a nonlinear model of tho F-GC 
aiccroft, there are sovernl kinds of probabilltloe 
responacs that one can guess boforo simulation, 
such os during tho transient period the MMAC al- 
gorithm tends to pick F/C's which aro mismatch 
stable while during the equilibrium condition the 
8* doi.ilnant behavior tends to occur, ilowovor, a 
rigorous statement on the precise nature of the 
probability rcBponnco is still an open question. 
Other thon tho model identification problem, the 
overall performance of the MMAC method system is 
very good. Although poorly selected F/C's in tho 
bank of Kalman filters may degrade the porformanciv 
tho MMAC method seems to stobillzo the oystem 
quite well. Tlic 8* dominant problem is basically 
tho same problem as in most parameter identifi- 
cation problems when there le lock of information. 
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Darling ratio for cloaod-Ioop short |M,riod 
polos SI a function of luxlnun pitch rate 
penalty, q . In (2.6) or (2.9) 


O..S 

Daitplng ratio 
for all sulisenlc 
comlltlons 

Danpirn ratio 
for all aupbrsenic 
coadltlnriV 

10,/V^ 

0.40B 

0,961 

B,/V^ 

0.590 

0.402 

C/V^ 

0,553 

0.440 

4»/V„ 

0.507 0.49D 

tabu: in 
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ri«ur* 5 ProUblUtlf* Veriui Timo (tru« r/C 09 
Ineludad) /rtm LongUudirUl pyn*nlcB 



rifiure 3 rrabAbllltlot Varsua TinsQ (tnm r/C 119 
net Included) from Lateral Dyniniea 



rlfuro 4 Uterel Acceleration vortun Tine Usinq 
Lateral infontvatlon only 



rlquta i Probabllltiei Voreue Tice (true r/c *19 
net included) froa Len9itudlnal Dynamic* 



rigura 7 Angle of Attack Varau* Tica) 
Longitudinal Infornition Only 




Appondlx A 


Wind Dloturbance Modal 


Ao rciurkod in soctlono 2 and 3^ a contlnuouB 
tlnu wind dlotutbanca inodal was Included for both 
the lateral and lon 9 ltudiiial dynamics « corro8-> 
ponding to a state variable w(t). In this appon- 
dlx wo give the mathematical details of this model 
which was kindly provided by Hr, J. Elliott of 
NASA/liRC, as 0 reasonable approximation to the von 
Karman model and the Haines approximation. It is 
Important to realize that the wind disturbance 
model changoB from flight condition to flight 
condition. The power spectral density of the 
wind disturbance is given by 


(A.l) 



( 


4 + 



) 


where h, the scale length, is 

I 200 ft at sea level 

(A, 2) L » . 2500 ft when altitude > 2500 ft 
I linoorly interpolated in between 

Vq is the speed of aircraft in ft/scc, U in 
rad/s oc, and 


I 6 ft/soc normal 

(A. 3) 0 X ^ 15 ft/scc in cumulus clouds 

I 30 ft/scc in thunderstorms 

To obtain a state variable model, a normal- 
ized statu variable w(t) (in rad) is used as the 
wind statu for both lateral and longitudinal 
dynamics. The state variable w(t) is the output 
of a first order system driven by continuous 
white noise 5(t-,) with zero mean. Thus the dynam- 
ics of the wind disturbance model are give by 


(A. 4) w(t) « -2(-^)w(t) + ■ ' C(t) 

f VroiT 

O 


where C(t} is zero mean white noise with unity 
covariance function 

(A. 5) E ^C(t)e{t>| - 6{t-T) 

The design was obtained for the intermediate case 
a X 15 (cumulus clouds) , 


For the longitudinal dynamics the wind state 
w(t) influences the dynamics in the same manner as 
the angle of attack. Thus, In the longitudinal 
state equations the Wind state w(t) enters the 
equations as follows 

I q(t) » + •,,''(t) 

(A. 6) ! v(t) a w(t) 

I . il u. 

( a(t) - ....... + a W(t) 


where a , a , a can be found frosi the open 
II II II 

loop longitudinal A aatrix 121. 

In the lateral dynamics the wind state w(t) 
Influences the dynamics in the same manner as thu 
sideslip angle. Thus, in the lateral state 
equations the wind state w(t) enters ths equa- 
tions as follows 

I p(t) - ♦ e|,w(t) 

(A. 7) V r(t) - + 

, I §(t) - 

whore a , a , a can be found from the open 

II Ji II 

loop lateral A matrix (2] . 


